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ABSTRACT 

We study the structural evolution of massive galaxies by linking progenitors and descendants at a 
constant cumulative number density of n c = 1.4 x 10~ 4 Mpc~ 3 to z ~ 3. Structural parameters were 
measured by fitting Sersic profiles to high resolution CANDELS HST WFC3 J125 and -ffieo imaging 
in the UKIDSS-UDS at 1 < z < 3 and ACS I S u imaging in COSMOS at 0.25 < z < 1. At a given 
rcdshift, we selected the HST band that most closely samples a common rest-frame wavelength so as 
to minimize systematics from color gradients in galaxies. At fixed n c , galaxies grow in stellar mass 
by a factor of ~ 3 from z ~ 3 to z ~ 0. The size evolution is complex: galaxies appear roughly 
constant in size from z ~ 3 to z ~ 2 and then grow rapidly to lower redshifts. The evolution in the 
surface mass density profiles indicates that most of the mass at r < 2 kpc was in place by z ~ 2, and 
that most of the new mass growth occurred at larger radii. This inside-out mass growth is therefore 
responsible for the larger sizes and higher Sersic indices of the descendants towards low redshift. At 
z < 2, the effective radius evolves with the stellar mass as r e oc M 2 ' , consistent with scenarios that 
find dissipationless minor mergers to be a key driver of size evolution. The progenitors at z ~ 3 were 
likely star forming thin disks with r e ~ 2 kpc, based on their low Sersic index of n ~ 1, low median 
axis ratio of b/a ~ 0.52, and typical location in the star forming region of the U— V vs. V — J diagram. 
By z ~ 1.5, many of these star forming disks disappeared, giving rise to compact quiescent galaxies. 
Towards lower redshifts, these galaxies continued to assemble mass at larger radii and became the 
local ellipticals that dominate the high mass end of the mass function at the present epoch. 
Subject headings: galaxies: structure — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

Massive galaxies in the nearby universe are gener- 
ally comprised of quiescent elliptical and SO galaxies. 
The formation of such massive galaxies has been an ac- 
tive area of study. Recent observations suggest that 
the properties of these quiescent galaxies (QGs) were 
much different at earlier times. For example, the sizes 
of QGs have been found to be muc h smaller, at fixed 
stellar mass, at h igh redshif t (e.g ., IDaddi et al.l 120051: 
Zirm et al.l [20071: iToft et all [20071: Ivan Dokkum et all 



2008; Cimatti et al. 2008 ; Ivan der Wei et al.l 2005 , 
Franx et al.l 120081: iWilliams I et al.l 120101: iNewman et al l 
20121) . This implies much higher mass densities within 

the effective radius for QGs at high redsh ift. T he size 

measurements are robust (Sz omoru et al.l 120101 120121) 
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and the stellar mass measurements are in good agree- 
ment with dynamical mass est i mates (|Cenarro fc Truiillol 
20091: iCappellari et all 120091: Ivan Dokkum et al.l 120091 



van de Sande et al.l 120111 ) ."confirming the dense nature 



of QGs at high redshift. Several mechanisms have been 
prop osed to explain the g rowth in sizes among QGs (see, 
e.g.. IHopkins et aLll2010D . Recent discussions have cen- 
tered on the relative importance of major and minor dis- 
sipationless mergers with the latte r favored to be the pri- 
mary channel for size growth (e.g. JBournaud et al.ll2007t 
INaab et al.ll200l iBezanson et al.ll2009b IHilz et al.ll2012l ). 
While evidence is emerging that th ese compact QGs rep- 
resent the cores of local ellipticals (IBezanson et al1l2009t 
IHopkins et "all 120091: Ivan Dokkum et all 120101) . the pro- 
genitors of these compact QGs at even higher redshifts 
remains a mystery. 

Selecting galaxy samples at or above a fixed stellar 
mass limit has provided important insight into the evo- 
lution in properties for such populations. However, the 
connection to the evolution of a typical galaxy over cos- 
mic time is not straightforward given that galaxies grow 
in stellar mass due to in situ star formation and merg- 
ing: the progenitors of galaxies that lie just above a stel- 
lar mass limit at low redshift would not be counted in 
a census of high redshift galaxies above the same mass 
limit since they were likely to be less massive. Mea- 
suring the structural evolution of a galaxy as it grows 
in mass therefore requires a method for linking its pro- 
genitors and descendants over cosmic time. One such 
method involves selecting galaxies at a constant cumu- 
lative number density. The basic principle behind this 
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method is that the rank ordering of galaxy masses does 
not change drastically over time. Therefore, if one selects 
the 10th most massive galaxy in a comoving volume at 
z ~ 3, it is still likely to be approximately the 10th most 
massive galaxy in that comoving volume at z ~ 0, but 
with a higher overall mass due to star formation and 
merging. This technique presents a complementary ap- 
proach to mass-selected studies. A number density se- 
lection has been used in other recent works to study the 
structural properties (Ivan Dokkum et al1l2010f) . star for- 
mation histories (SFHs. [ Papovich et al.ll2011l ) , and mass 
growth (iBrammer et alJ 120111) of gal axies over cosmic 
time (see also. lLoeb fc Peebles! 120031) . The study by 
van Dokkum et all (|2010D was carried out with ground 
ba sed imaging. In this wo rk, we build on the results 
in Ivan Dokkum et all (|2010| ) and utilize high resolution 
Hubble Space Telescope (HST) imaging, allowing us to 
measure structural properties more accurately and to 
push to z ~ 3. At these high redshifts, we also iden- 
tify the progenitors of ~ 2M* galaxies, which in the lo- 
cal universe are bulge dominated, quiescent systems with 
large effective radii. 

The layout of this paper is as follows. In Section [2] 
we discuss the data used for the study. In Section [3] we 
discuss the relevant derived quantities. In Section U we 
examine the structural assembly of a galaxy as it grows 
in time to become a ~ 2M* galaxy by z ~ 0. Our re- 
sults are further discussed in Section[5]and we summarize 
our findings in Section |6l For completeness, we briefly 
overview the structural properties of QGs and star form- 
ing galaxies (SFGs) above a constant stellar mass limit in 
the Appendix, as this provides an alternative view to the 
number density selection in the main part of the paper. 

We assume a cosmology with Hq = 70 km s _1 Mpc -1 , 
Qm = 0.30, and O a = 0-70- Stell ar masses are based 
on a Chabrier IMF (jChabrierll200l. Adopting IMFs of 
differ ent forms below 1 Mq (e.g. . Ivan Dokkum fc Conrovl 
I2010D would lead, to first order, to an overall scaling of 
the stellar masses and otherwise identical results. All 
magnitudes are given in the AB system. 

2. DATA 

2.1. UDS-CANDELS: 1< z < 3 

We study the structural properties of massive galax- 
ies at high redshift (1 < z < 3) using a combina- 
tion of data f rom the UKIRT Infr ared Deep Sky Sur- 
vey (UKIDSS. lLawrence et al.ll2007f ) and the Cosmic As- 
sembly Near- infrared Deep Extra galactic Legacy Survey 
(CANDELS, iGrogin et al1l20ll . We employ a multi- 
wavelength dataset (u*BVRi'z JHK and Spitzer IRAC 
3.6 /Ltm and 4.5 /zm) in the Ultra-Deep Survey (UDS) 
field of UKIDSS. This is one of the deepest wide-field 
near-IR datasets available, making it ideal for construct- 
ing samples of galaxies to relatively low stellar masses 
at high redshift. The obs ervations and d a ta reduction 
are de scribed in de t ail in iWilliams etldl (|2009l 12010ft 
and in lQuadri et al.l ()2012l ) and briefly summarized here. 
The UKIDSS UDS DR8 data include JHK, which have 
5cr limiting depths in D — l'/8 apertures of 24.9, 24.1, 
and 24.5 AB mag, respectively. The BVRi'z' imaging 
were obt ained as part of the Subaru-AM M Deep Survey 
(SXDS, ISekieuchi fc SXDS Team! l200l . The u* data 
was obtained with MEGACAM on CFHT (PI: O. Al- 



maini). The Spitzer IRAC 3.6 fim and 4.5 fim data were 
obtained as part of the Spitzer-XJDS Survey (SpUDS, PI: 
J. Dunlop). An updated catalog will be presented in 
detail in Williams et al. (2012, in preparation). Ob- 
jects were detected in the if-band and photometry car- 
ried out in the other bands with matched apertures. The 
IRAC photometr y was measur e d usin g the PSF convolu- 
tion procedure of lLabbe etall ((2006). The IRAC bands 
in the UDS are crucial for determining photometric red- 
shifts, stellar masses, and rest-frame optical/near- IR col- 
ors of galaxies at high redshift. An analysis of simulated 
number counts indicates that the catalog is > 90% com- 
plete at K tot = 24.0 AB mag. 

To measure structural parameters we utilize HST 
imaging in the UDS that was acquired as part of CAN- 
DELS^t fGrogin et al.l l2 011ft. The observations and data 
reduction are presented in iKoekemoer et al.l (|2011h . We 
employ ul.O of the publicly available WFC3 J125 and 
iJi6o mosaics. The PSF FWHM in these two bands are 
0'.'12 and 0'.'18, respectively. The CANDELS HSTWFC3 
data cover roughly ~ 0.06 deg 2 of the ~ 0.65 deg 2 UDS 
field with multi-wavelength coverage. 

2.2. COSMOS- ACS: 0.25 < z < 1 

In order to assemble a large sample of massive galax- 
ies at low redsh ift (0.25 < z < 1), w e use data in the 
COSMOS field (IScoville et alJ l2007blah . This ~ 2 deg 2 
field benefits from multi-wavelength imaging spanning 
the UV to IR, i ncludi ng wide-field HST ACS cover- 
age. Illbert etaLl fj2009f ) assembled a photometric cata- 
log in COSMOS incorporating CFHT u* and K, Sub- 
aru BVgriz as well as 12 intermediate optical bands 
(IA427, IA464, IA484, IA505, IA527, IA574, IA624, 
IA679, IA709, IA738, IA767, IA827), UKIRT J, all four 
IRAC channels, and GALEX NUV and F UV. Some of 
these d ata products are fu rther described in lCapak et alJ 
(2007). Illbert et al.1 (|2010t l used the photometry in COS- 
MOS to derive stellar mass estimates, which could in 
principle be used in our analysis. However, for con- 
sistency with the data preparation and stellar masses 
computed in the UDS (see Section [3j), we use a K s - 
selected catalog, which incorporates most of the data in 
COSMOS described above but reconstructed in a man- 
ner consistent with what was done in t he UDS (e.g., 
IWilliams et "all 12009. lOuadri et all |2012| ). A detailed 
account of the observations, data reduction, and cat- 
alog construction will be described in a future paper 
(Muzzin et al., in preparation). We briefly remark on 
the most relevant points of the catalog here. Objects 
were detected in the if s -band, which reached a 5a detec- 
tion limit of 23.85 AB mag for a D — 2" aperture. An 
analysis of simulated number counts indicates that the 
90% completeness is K s ^ ot = 23.5 AB mag. 

The primary purpose of using data in the COSMOS 
field is the wide-field HST ACS igi4 imaging, which we 
use to measure structural parameters (see Section [3]). 
In this work, we utilize the v2.0 ACS J si4 imaging 
(|Koekemoer et al.l 120071 : iMassev et alJ 12010ft . The typ- 
ical PSF FWHM is O'.'l. The ACS imaging employed 
here covers roughly ~ 1.3 sq deg of the COSMOS field. 



http: / /candcls. ucolick.org/data_access/Latest_Release. html 
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3. ANALYSIS 

3.1. Photometric Redshifts, Stellar Masses, and 
Rest-frame Colors 

While the datasets used in this work are assembled 
from two different fields, the analysis carried out on the 
data is uniform. Photometric redshifts in both COSMOS 
and UPS were measured with EAZY (jBrammer et al.l 
2008). In the UDS, a comparison to spectroscopic red- 
shift measurements suggests that the photometric red- 
shift uncertainties are c 2 /(l + z) ~ 0.022 at 1 < z < 1.5. 
At higher redshifts, where the majority of our sam- 
ple in the UDS lies, iQuadri et all ((2012) find that the 
photometric redshift uncertainties are larger based on 
differences in the photom etric redshifts of close pairs 
(IQuadri fc Williamsl[20Tol) . In COSMOS, a similar com- 
parison to spectroscopic redshifts at 0.25 < z < 1 in- 
dicates that the photometric redshift uncertainties are 
<J Z /(1 + z) ~ 0.01. Rest-frame U — V and V — J colors 
were also computed with EAZY. 

Stellar masses were computed with the SED fit- 
ting code FAST dKriek et all [2 009) usi ng ex ponentially 
declining SFHs. iBruzual fc Chariot! (|2003l . hereafter, 
BC03) stellar population synthesis (SPS) models with 
a Ch abrier IMF w e re use d in the SED fitting. In the 
UDS, IQuadri et al.l (|2012f ) parameterize the stellar mass 
limit as a function of redshift as M\i m — 9.4 + 1.21n(z). 
Using the same tec hnique to compute the mass limit as 
IQuadri et al.l (|2012f) , but focusing solely on galaxies near 
z = 3, we find that a mass limit of M ~ 10 la6 Mq en- 
compasses ~ 95% of galaxies. We adopt this value as our 
stellar mass limit in the UDS at z = 3. As will be made 
clear in Section SI galaxies in our COSMOS sample at 
0.25 < z < 1 are well above the limiting stellar mass at 
z = 1 of M~ 10 9 - 4 Mq. 

3.2. Structural Parameters 

Structural para meters were obtained with GALFIT 
(|Peng et al.l 120021) . which provided measurements of 
Sersic indices (n), effective radii (r e ), and axis ratios 
(b/a). The effective radii reported here are circularized, 
r e = y/ab. In order to carry out these measurements at 
the same rest-frame wavelengths, we use the HST imag- 
ing that is closest in rest wavelength to Ao = 5160 A. 
This choice of Ao represents a tradeoff between (i) the de- 
sire to push to high redshifts, and (ii) the need to probe 
rest-frame light as redward as possible from the 4000 A 
break, beyond which the SEDs of galaxies are generally 
smooth and serve as a closer tracer of the stellar mass. 
As a result, in COSMOS, we use ACS igi4 imaging at 
0.25 < z < 1, while in the UDS we use WFC3 J 125 
imaging at 1 < z < 1.76 and WFC3 H\eo imaging at 
1.76 < z < 3. For a given object, the nearest star is cho- 
sen to serve as the PSF model when running GALFIT. 
Nearby objects were masked. No constraints were placed 
on the range of Sersic indices. The axis ratios were con- 
strained to have 0.1 < b/a < 1. The semi-major axis, a, 
was constrained to be smaller than the box size. 

Our results and conclusions are dependent on the qual- 
ity of the structural parameters measured from GAL- 
FIT. We ran simulations to test the reliability of these 
measurements at the highest redshifts, where objects be- 
come faint in the WFC3 Hi 60 imaging. We created 



10, 000 mock galaxies from a range of Sersic models, vary- 
ing r e , n, b/a, and the magnitude. These models were 
added to different regions of the WFC3 Hi 60 imaging 
that contained blank sky. We then processed the mock 
images in the same manner as used to derive our mea- 
surements above. The sample in our highest redshift bin 
(2.5 < z < 3), defined in Section HI reaches magnitudes 
of ifieo ~ 23.6 AB mag and 24.3 AB mag at the 50th 
and 90th percentile, respectively. Our simulations show 
that we recover r e , n, and b/a to precisions of ~ 11%, 
5%, and 2% for H wo ~ 23.6 AB mag and - 23%, 12%, 
and 5% for -Hi6o ~ 24.3 AB mag. Systematic offsets are 
< 1% for all three parameters. Effective radii as small 
as 0'/06 (i.e., less than the FWHM/2), or r e ~ 0.5 kpc at 
2 < z < 3, are recovered to similar prec isions as noted 
abov e, consistent with other works (e.g.. INewman et al.l 
120121) . Our results and conclusions based on structural 
parameters measured from GALFIT are therefore not 
strongly biased at high redshift in any way. 

As one might expect at higher redshifts, the higher res- 
olution HST imaging reveals that a small portion of the 
objects detected in the ground based if-band imaging are 
comprised of two or more objects blended together. We 
remove galaxies in the UDS from our analysis (< 10%) 
if nearby objects contribute more than 10% of the total 
flux within the D = l'/8 color aperture, as the photomet- 
ric redshifts, rest-frame colors, and stellar masses are not 
as reliable. We note that this procedure may introduce 
a small bias against galaxies in close pairs. 

4. THE ASSEMBLY OF MASSIVE GALAXIES 

We define our sample in this section by selecting 
galaxies at a fixed cumulative number density, n c . We 
then study the structural properties of galaxies in nar- 
row mass bins at this value of n c at different redshifts. 
The analysis that follow s is overall similar to that of 
Ivan Dokkum et all ()2010D but differs in a few aspects. 
First, we study galaxies at a constant cumulative number 
density. Second, our analysis extends to z ~ 3. Finally, 
our work employs HST imaging as opposed to ground 
based, allowing us to more accurately measure struc- 
tural properties, especially at high redshift. The deep 
HST imaging also allows us to characterize the proper- 
ties of individual galaxies as opposed to s tacks, which 
were employed in Ivan Dokkum et all Ip010[ ). 

4.1. Selection at a Constant Cumulative Number 
Density 

Figure [TJj shows the cumulative number density of 
galaxies at different redshifts. These curves were derive d 
using the mass functions of Marchcsi ni et al.l (|2009f> . 
which were computed with data from several fields and 
as a result minim ized uncertainties from cosmic vari- 
ance. In addition, Marchcs ini et al.l (|2009|) carefully ac- 
counted for completeness limits in stellar mass in de- 
riving their mass functions. We used best-fit Schechter 
parameters from Set 7 in that work, which used so- 
la r metallicity BC03 m odels with a Chabrier IMF and 
a ICalzetti et al.l ((2000) dust law to determine stellar 
masses, similar to the SFHs employed in this work. We 
derived mass functions for our redshift bins by interpo- 
lat ing between the best-fit S chechter functions computed 
bv iMarchesini et al.l (|2009fl . which cover 1.3 < z < 4. 
In carrying out the interpolation, we also included the 
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ICole et all (12001 m ass fu nction at z ~ 0.1, as reported 
in lMarchesini et al.l (|2009h but scaled to a Chabrier IMF. 
We then integrated these mass functions, $(M), to de- 
termine the cumulative number density of galaxies as a 
function of stellar mass at different redshifts: 



n{> M) 



$(M)dM 



(1) 
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Figure 1. (a) Cumulative number density of galaxies at a given stellar mass for different redshifts, derived from the mass functions of 
IMarchc sini et al. ( 2009|). At a fixed cumulative number density of n c = 1.4x 10 — 4 Mpc~ 3 (dashed black line) we determine the corresponding 
stellar mass for a given redshift bin (dotted vertical lines), (b) Stellar mass vs. redshift for galaxies selected at n c = 1.4 X 10 — 4 Mpc -3 . 
The solid curve represents a second order polynomial fit and is given by Equation [2] A galaxy with a stellar mass of M 5 X 10 10 Mq 
at 2 = 2.75 grows by a factor of ~ 3 in mass by z = 0.375. For a given redshift, we study the structural properties of galaxies at 
n c = 1.4 X 10 — 4 Mpc -3 by selecting objects in a narrow mass bin around the predicted stellar mass from Equation[2] 

mass function, in practice this results in selecting galax- 
ies at (logM„ c /M Q )t^ 5 . The bin size is broad enough 
to allow for robust measurements of median structural 
parameters. In order to avoid confusion, we emphasize 
that at a given redshift we are not selecting all galaxies 
with masses above the mass limit implied by the given 
value of n c , but instead, we are selecting galaxies in a 
narrow mass bin at the mass determined by n c . 

Equation [2] indicates that galaxies at n c = 1.4 x 
10 -4 Mpc -3 grow by a factor of ~ 3 from z = 2.75 to 
z ~ 0, resulting in a galaxy at low redshift with a stellar 
mass of M ~ 1.5 x 10 11 M Q (i.e., ~ 2M*). From z = 2 to 
z = 0.1, Equation [5] predicts that the stellar mass grows 
by a factor of ~ 2, wh i ch is y ery similar to what is found 
in Ivan Dokkum et al.l {2010). We note that the stellar 
mass growth inferred from Equation [5] is less than what 
is predicted from abunda nce matching techniques (e.g., 
iConrov fc Wechslerl 120091 ). though the latter analysis is 
quite uncertain at z > 1. At z ~ 0, where the obser- 
vations are more robust, abundance matching indicates 
that a ~ 2M* galaxy occupies a dark matter halo of mass 
M ~ 2 x 10 13 M(D, which is typica l of ga laxy groups. 

Bot h Ivan Dokku m et al.l (|2010D and IPapovich et al.l 
(|201lD show with simulations that in selecting galaxies 
at a fixed number density, the completeness fraction de- 
clines with cosmic time, meaning that some of the objects 
selected in a given number density bin at high redshift 
are no longer found in that bin at lower redshift. Con- 
taminants from other number density bins also enter the 
sample. However, most of the contaminants scatter in 
from neighboring bins and likely display properties that 
are very similar to those of galaxies in the number density 
bin of interest. 

Fi nally, we note that altho ugh we use mass functions 
from IMarchcsi ni et al.l (|2009f ). which are based on dif- 
ferent data from what is employed here, we arrive at 
qualitative and quantitative conclusions that are quite 



We chose a cumulative number density for our study of 
n c = 1.4 x 10~ 4 Mpc~ 3 (dashed line) as this value repre- 
sents the number density of galaxies with stellar masses 
slightly above the stellar mass limit at z = 3 (M ~ 
iqio.6 Mq). As seen in Figure [TJi, the selected value 
of n c intersects the cumulative number density curves at 
lower redshifts at higher stellar masses, tracing out the 
mass growth at that particular number density. The cor- 
responding stellar masses at n c = 1.4 x 10~ 4 Mpc -3 are 
shown as a function of redshift in Figure [2). In order 
to quantify the redshift dependence we fit a second order 
polynomial to these datapoints, resulting in the following 
relation between stellar mass and redshift: 



logM„ c /M = 11.19 - 0.068z - 0.040z : 



(2) 



The scatter about this relation is only a = 0.0046 dex, 
suggesting that this parametrization is adequate for the 
redshift range studied here. Note that this scatter does 
not reflect the systematic uncertainty in measuring stel- 
lar masses of galaxies at high redshift , which can be sub- 
stantial (e.g. . IMarchcsi ni et alj 12009). For a given red- 
shift, we study the properties of galaxies within a bin of 
size ^0.3 dex in stellar mass centered on the predicted 
mass from Equation O The actual boundaries of the bin 
are adjusted such that the median mass is close to the 
value given by Equation [2] Given the steepness of the 
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Figure 2. Rest-frame U — V vs. V — J color for galaxies at a fixed cumulative number density, n c = 1.4 X 10 — 4 Mpc - 3 , in COSMOS 
(0.25 < z < 1) and the UDS (1 < z < 3). Galaxies were selected in narrow mass bins at each redshift around the predicted stellar mass 
in Equation [2) Note that the number density selection results in a selection of more massive galaxies towards low redshift (e.g., Figure [TJ- 
The gray points indicate the parent sample of galaxies while objects color-coded by Sersic index indicate the sample with HST imaging 
(note that the Sersic indices were not bound to the range shown). The sample size of galaxies with measured structural parameters is 
indicated in the bottom right of each panel. The solid black line shows the division between QGs (top left) and SFGs (bottom right). At 
z > 2, the progenitors of ~ 2M* galaxies were star forming disks. By z ~ 1.5, many of these compact star forming disks have disappeared 
from the sample, while compact QGs have emerged. At z ~ 0.375, the assembled ~ 2M* galaxies resemble bulge dominated, quiescent 
systems with large effective radii. 

similar to what was found in Ivan Dokkum et al.l (2010). 
This suggests the systematic uncertainties as a result of 
this choice are minimal. 
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4.2. Star Formation Properties 

We first examine how the star formation properties 
of galaxies have evolved since z ~ 3. As is well 
known, galaxies can be classified in two distinct cate- 
gories: star forming an d quiescent, at least out to z ~ 3 
(Whitakc r et al.ll2011f) . Figure [5] shows the rest- frame 
U — V vs. V — J colors of galaxies in different redshift 
bins selected at n c = 1.4 X 10~ 4 Mpc" 3 . This UVJ dia- 
gra m is commonly used to separate QGs from SFGs (see, 
e.g.JLabbe et all2006tlWuvts et al.ll2007t lWilliams et al.l 
l2009HPatel et al.ll201lLl2012D . It ispreferred over a color- 
magnitude or a color-mass selection because of its ability 
to separate red galaxies that are quiescent from reddened 
SFGs. Shown in each redshift panel are galaxies within a 
narrow mass bin around the predicted stellar mass from 
Equation [5] for galaxies at n c . The typical mass is there- 
fore increasing towards low redshift. The subset of galax- 
ies with measured structural parameters from the HST 



ima ging are color-coded a ccording to their Sersic index. 
The I Williams et aLl (|2009T > boundary distinguishing QGs 
(top left) from SFGs (bottom right) is shown for each 
redshift bin. At z > 2, we use the boundary c ondition 
defined at 1 < z < 2 since IWilliams et al.l (|2009f ) provide 
selection criteria up to those redshifts. We also slightly 
modify the diagonal boundary at 0.5 < z < 1 to better 
accommodate the COSMOS data using the following re- 
lation: U-V> 1.08 x {V- J) +0.43. Based on the UVJ 
selection, Figure [3] shows the fraction of QGs and SFGs 
as a function of redshift. The error bars in Figure [3] at 
z > 1 are larger due to the smaller sample size of the 
UDS data. Below z < 1, the wide area COSMOS data 
allow us to better constrain the properties of the most 
massive galaxies. The QG fraction increases from ~ 23% 
at z ~ 2.75 to - 89% at z ~ 0.375. We note that our 
quiescent and star forming fractions are in good agree- 
ment with those of Bramm er et al.l (|2011l) who reported 
on galaxies to z ~ 2 (open symbols in Figure [3]) . 

Together, Figures [2] and [3] show that the progenitors 
of nearby massive galaxies at z ~ 3 were likely to be 
star forming given that their colors coincide with those 
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Figure 3. Fraction of UVJ classified quiescent galaxies (QGs, 
solid red circles) and star forming galaxies (solid blue squares) vs. 
redshift for galaxies with measured structural parameters selected 
at a constant cumulative number density of n c = 1.4x10 Mpc -3 . 
The star forming fraction is simply the complement of the quies- 
cent fraction. The la error bars are computed assuming a bino- 
mial distribution. The open symbol s represent the values for the 
appropriate mass and redshift from Brammcr ct al. (2011). The 
change in the proportion of QGs towards low redshift for galaxies 
at n c is dramatic, increasing from ~ 23% at z ~ 2.75 to ~ 89% at 
2 ~ 0.375. At z ~ 3, most of the progenitors of massive galaxies 
were star forming. 

of SFGs in UVJ color space. From z ~ 3 to z ~ 1.5 
there appears to have been a substantial buildup in the 
population of QGs, which become the dominant popula- 
tion at z < 1.5 for galaxies selected at n c . At z < 1.5, 
the population of QGs continues to grow. We examine 
how the structural properties have evolved in the next 
section. 

4.3. Structural Evolution 

We now turn to the structural properties of massive 
galaxies as they grow in time. The general proper- 
ties of galaxies selected at n c are summarized in Ta- 
ble Q] Figure 0] shows the evolution of the median ef- 
fective radius, Sersic index, and axis ratio for galaxies 
at n c = 1.4 X 1CP 4 Mpc -3 . We include in these figures 
a sample of SDSS galaxies at z = 0.06 from Szomoru 
et al. (2012, in preparation) selected at the appropriate 
mass from Equation [2] The SDSS data is not used in 
any of the fits that follow. Below z < 2, where QGs be- 
come the dominant population in the sample, Figure |4ji 
shows that the effective radius increases substantially 
from r e ~ 2 kpc at z ~ 2 to r e ~ 7 kpc at z ~ 0. 
The size evolution at 0.25 < z < 2 follows 



= (9.3 ±1.0) kpcx (1 + 2) 



-1.1±0.2 



(3) 



wi th the exponent being ver y similar to what was found 
in Ivan Dokkum et all (|201Ct ) over roughly the same red- 
shift range. A striking feature in Figure [4h is the lack of 
evolution at 1.5 < z < 3 in the median effective radius. 
We investigate this further by showing the evolution of 
QGs and SFGs separately in Figure |4j The constant 
median r e arises because SFGs are larger than QGs and 



their relative abundance changes as a function of redshift. 
Above z > 3, the size evolution is likely determined al- 
most solely by SFGs since they become an overwhelming 
majority of the population. We theref ore expect the sizes 
of galaxies to decr ease above z > 3 (jOesch et alJl20ld 
iMosleh et alJl2012l see, e.g.,) for samples selected at n c . 
We can test whether the apparent constant value of r e at 
1.5 < z < 3 is a generic feature or a consequence of the 
particular value of n c selected for our study. At lower 
values of n c (i.e., higher masses at a given redshift), we 
find that r e can increase gradually from z ~ 3 to z ~ 2. 

The Sersic index determines the distribution of light 
and hints at the presence of a bulge or disk. The median 
Sersic index in Figure |4f) increases from n ~ 1 at z = 
2.75 to n ~ 6 at z ~ 0. The Sersic index evolution at 
0.25 < z < 3 can be characterized by 



(6.7+0.5) x (1 + z) 



-0.9±0.1 



(4) 



which is again very sim ilar to what was found in 
van Dokk um et al.l (|2010f ). The Sersic index evolution 
indicates that while most of the stars in ~ 2M* galaxies 
in the nearby universe are distributed in a bulge, the stars 
in their progenitor galaxies at z > 2 were distributed in 
structures resembl ing exponential disks. We note that 
iWuvts et all (j201lD also find that SFGs at high redshift, 
which represent the majority at n c , generally have low 
Sersic indices around n ~ 1. The fact that QGs at 
2 < z < 3 in Figure [2] g enerally have hig her Sersic in- 
dices than SFGs (see also. lBell et al.l [2012). as is also the 
case at lower redshifts, further suggests that our Sersic 
profile fitting measurements are not significantly biased 
by the limiting depth of the HST imaging for higher red- 
shift galaxies. In the Appendix, we show this to also be 
the case with a much larger, stellar mass limited sample 
(Figure EJ. 

While the Sersic indices can be suggestive of a bulge 
or disk component, the axis ratio distribution provides 
a better constraint on the shapes of galaxies. Owing to 
the high resolution of the HST imaging, we can examine 
the axis ratios of galaxies selected at n c to z ~ 3. The 
median axis ratio of galaxies at n c = 1.4 x 10~ 4 Mpc -3 
has increased significantly since z ~ 3. At z = 2.75, the 
typical axis ratio is b/a ~ 0.52, a low value that is in- 
dicative of a distribution of randomly oriented thin disks. 
Meanwhile, at z — 0.06, the axis ratio is b/a ~ 0.76, 
closer to what is expected for elliptical galaxies. This 
value is in good ag reement with SPSS studies of massive 
QGs at z = 0.06 ijvan der Wei et al.ll2009t iHolden et al.l 
120121 ) . At the highest redshifts (2.5 < z < 3), the 
residuals to the single component Sersic profile fits are 
smooth and visual inspection of these residuals suggests 
that the lower axis ratios at high redshift are not gener- 
ally driven by multiple components. In addition, vary- 
ing the number density selection to lower values, such as 
n = 10 -4 Mpc~ 3 , does not impact the general decreasing 
trend of the median axis ratios towards higher redshifts. 
Finally, as seen in Figure [2b the axis ratios of SFGs are 
generally lower than that of QGs at a given redshift and 
SFGs increasingly become the dominant population at 
z > 1.5. The axis ratios therefore also indicate, in addi- 
tion to the Sersic indices, that the stars in the progenitors 
of ~ 2 M* galaxies were distributed in disks at z ~ 3. 

Figure [5] shows example postage stamps of galaxies se- 
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Figure 4. Evolution in the structural properties of galaxies selected at a constant cumulative number density of n c = 1.4 X 10 — 4 Mpc — 3 . 
(a) Effective radius vs. redshift. The shaded region indicates where r e < FWHM/2. (ft) Sersic index vs. redshift. (c) Axis ratio vs. 
redshift. The black circles represent the median for the full n c selected sample while the red and blue data points represent the median 
values for the QG and SFG sub-populations, respectively (offset in redshift for clarity). The gray filled circles at z = 0.06 represent the 
median values for an SDSS sample at z = 0.06 from Szomoru et al. (2012, in preparation). The typical size of a galaxy at n c increases by 
a factor of ~ 3 — 4 since z ~ 3, with most of this change occurring at z < 2. The apparent constant size at 1.5 < z < 3 is a consequence 
of the changing mix of QGs and SFGs combined with the overall growth of galaxy masses and sizes over this redshift range. The Sersic 
indices increase from n ~ 1 at high redshift to n ~ 6 at low redshift. This suggests that most of the stars in galaxies at n c were distributed 
in a disk at z ~ 3, while at low redshift they are distributed in a bulge. The increasing axis ratios towards low redshift further support 
this view. The median axis ratio of b/a ~ 0.52 at z ~ 2.75 is close to what is expected for randomly oriented thin disks, while at z ~ the 
axis ratios are more indicative of spheroidal systems. 




Figure 5. Example postage stamps for galaxies selected at a constant cumulative number density of n c = 1.4 X 10 — 4 Mpc — 3 at different 
redshifts. At a given redshift, objects were selected to have properties that follow the general trends seen in Figure|4] Each postage stamp 
is 30 kpc on a side and is rotated such that the major axis is aligned horizontally. The PSF FWHM is indicated by the circle in the bottom 
right of each panel and the effective radius (in kpc), Sersic index, and axis ratio (q) are given in the bottom left. The half-light ellipse is 
shown in white. Towards lower redshift, more light is added to the outer parts, leading to the larger sizes, Sersic indices, and the buildup 
of stellar mass. Towards higher redshifts, the median axis ratio declines, suggestive of randomly oriented disks. 



lected at n c at different redshifts. For illustrative pur- 
poses, we selected galaxies with structural parameters 
that follow the general trends seen in Figure [4] but with 
more consideration for the trend in axis ratios. Each 
postage stamp is 30 kpc on a side and note that the size 
of the galaxy in each redshift bin is larger than the PSF. 
The relative sizes between redshift bins are therefore easy 
to compare in this figure when paired with the indicated 
half-light ellipses. Towards lower redshifts, especially at 
z < 2, Figure [5] shows how light is added to the outer 
parts, leading to the increasing size and Sersic index of 
galaxies selected at n c . At higher redshifts, the axis ra- 
tios decrease, suggesting a larger contribution from disks. 

Synthesizing the structural information above with the 
star formation activity discussed in Section 14.21 we see 
that the progenitors of ~ 2M* galaxies at z ~ 3 were star 
forming disks with r e ~ 2 kpc. By z ~ 1.5, many of these 
star forming disks disappear and give rise to a population 
of compact QGs. By z ~ 0, these compact QGs evolved 
into the large-sized, bulge dominated, quiescent ~ 2M* 
galaxies. 



4.4. Mass Assembly 

A more intuitive view of how the change in structural 
properties has impacted the evolution of galaxies at n c — 
1.4 x 1(T 4 Mpc" 3 is given by the evolution of the mass 
surface density profiles in Figure |BJi. In order to compute 
these profiles, we first used the best fitting Sersic index 
and effective radius for each galaxy to determine the light 
profile using the standard formula for a Sersic profile 

E(r) = £ e exp(-& B [(r/r e ) 1/n -l]) (5) 

where S(r) is the surface brightness at radius r, r e 
the half-light radius, E e the surface brightness at r e , 
n t he Sersic index , and b n a constant that depends on 
n. ISzomoru et al.l (|2012[ ) find that the surface bright- 
ness profiles of galaxies at high redshift generally follow 
Sersic profiles quite well. For simplicity, we therefore use 
the analytic representation of the profile for each galaxy. 
These light profiles were converted into mass profiles by 
normalizing the integrated light in the Sersic profile to 
the stellar mass of each galaxy. Note that this conversion 
neglects radial gradients in the mass-to-light ratio. For 
a given redshift bin, the median of the mass profiles was 



Table 1 

Properties of Galaxies Selected at a Constant Cumulative Number Density of 
n c = 1.4 X 1(T 4 Mpc~ 3 



Redshift 




Mass b 


Quiescent c 




Sersic c 


b/a f 


Range 




log M/Mq 


Fraction 


(k P c) 


Index 




0.25 < z < 0.5 


120 


11.16 


0.89 ±0.03 


6.1 ±0.3 


4.6 ±0.2 


0.74 ±0.01 


0.5 < z < 1 


118 


11.12 


0.77 ±0.02 


5.1 ±0.2 


4.1 ±0.09 


0.71 ±0.01 


1< z < 1.5 


20 


11.05 


0.65 ± 0.1 


3.2 ±0.5 


3.2 ±0.5 


0.64 ± 0.05 


1.5 < z < 2 


11 


10.95 


0.44 ± 0.08 


2.4 ±0.3 


2.5 ±0.3 


0.57 ±0.07 


2 < x < 2.5 


37 


10.84 


0.32 ±0.08 


2.3 ±0.3 


1.8 ±0.4 


0.68 ±0.06 


2.5 < z < 3 


35 


10.70 


0.23 ±0.07 


2.3 ±0.3 


1.5 ±0.3 


0.52 ±0.03 



a Number of galaxies in the sample at the given redshift. 

b Stellar mass of galaxies at n c = 1.4x 10~ 4 Mpc -3 for a given redshift (see Equation[2]l • 
c Fraction of galaxies that are quiescent based on UVJ selection. 
d Median effective radius for galaxies at n c . 
e Median Sersic index for galaxies at n c . 
f Median axis ratio for galaxies at n c . 




Figure 6. (a) Stellar mass surface density profiles of galaxies selected at a constant cumulative number density of n c = 1.4 X 10 -4 Mpc -3 
for different rcdshifts. At a given redshift, light profiles were derived for each galaxy in the 0.3 dex mass bin based on the best fitting Sersic 
index and effective radius. These light profiles were then normalized to the stellar mass of each galaxy and then median combined. The 
dotted portion of each profile indicates where the bootstrapped uncertainty of the median is greater than 20%. The gray shaded region 
extends to the radius that corresponds to the maximum PSF FWHM/2 for the full sample (occurs at z = 1.76). (b) Cumulative stellar mass 
at a given radius relative to the total mass within r < 100 kpc for a galaxy at z ~ 0. The mass profiles overlap at small radii suggesting 
very little mass growth in the inner parts of a galaxy at n c , while at larger radii there appears to be a substantial buildup of mass with 
cosmic time. 



computed at each radius resulting in the profiles shown 
in Figure Oi. The uncertainty in the median of the mass 
profile at a given radius was computed by bootstrapping 
the sample. 

A naive interpretation of Figure [H> would be that the 
bulges of galaxies grow in time given the increase in the 
Sersi c index, a crude proxy for the bulge-to-disk ratio 
(e.g. JLackner fc Gunnll2012l ). However, the mass profiles 
in Figure [6jz generally overlap at small radii and diverge 
at large radii, suggesting a buildup of mass in the outer 
parts of the galaxy with time. Figure© shows the cumu- 
lative proportion of mass assembled at different radii rel- 
ative to the total mass within r < 100 kpc of the median 
galaxy at z ~ 0.06. Roughly ~ 50% of the total mass of 
the galaxy is assembled within r < 7 kpc at z = 0.06, 
as expected given that r e ~ 7 kpc at that redshift. At 
z ~ 2.25, the assembled mass within r < 7 kpc is ~ 40% 



of the total mass at z = 0.06 indicating that much of the 
mass within r < 7 kpc was already in place ~ 10 Gyr ago. 
Note that the small sample in the 1 < z < 1.5 bin likely 
leads to this curve falling slightly above the < z < 1 
data at r < 10 kpc. 

In Figure we compare the mass growth between the 
central and outer regions of galaxies selected at n c = 
1.4 x 10~ 4 Mpc -3 . The total stellar mass as a function 
of redshift is shown by the black line, while the projected 
mass inside and outside of r — 2 kpc is given by the red 
and blue lines, respectively. These values are determined 
by integrating Equation [5] as follows: 

S(r)27rrdr (6) 
where r- ln and r out are the inner and outer radii enclos- 
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Figure 7. Projected stellar mass for different radial regions of 
galaxies selected at a constant cumulative number density of n c = 
1.4 X 10 — 4 Mpc - 3 . At small radii (r < 2 kpc, red line), most of 
the stellar mass was in place by z ~ 2. At larger radii (r > 2 kpc, 
blue line), there has been a substantial buildup of mass, fueling the 
overall mass growth of the galaxy. 
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Figure 8. Radius enclosing a fixed stellar mass as a function 
of redshift for a galaxy at n c = 1.4 X 10 — 4 Mpc -3 . Different 
colored lines indicate the evolution in the radius that encloses the 
given mass, ranging from 2 X 10 10 - 10 11 Mq. The lines therefore 
represent horizontal cuts in Figure |EJ) (ignoring a normalization 
factor). The lines of constant mass terminate in high redshift bins 
where the given mass had not yet been assembled. For reference, 
the effective radius for galaxies at fixed n c is indicated by the gray 
line. Below z < 2, the radius enclosing a given mass remains 
roughly constant, indicating that new stellar mass growth towards 
low redshift occurs at larger radii. From z ~ 3 to z ~ 2, the radius 
enclosing a given mass decreases but this trend may not be robust. 

ing the mass, M. For the central regions (n n = kpc, 
T'out = 2 kpc), the stellar mass appears to grow from z ~ 
3 to z ~ 2 but then levels off around - lfj 10 - 5 - 10 - 6 M . 
In contrast, in the outer regions (r; n = 2 kpc, r out = 
100 kpc) mass continues to build up over the entire red- 
shift range studied, growing by a factor of ~ 3. The 



stellar mass that has been added to the outer parts of 
galaxies over time is therefore the dominant source of as- 
sembled mass, as the central parts appear to have been 
assembled by z ~ 2. These results are in qual itative 
agreement with those of Ivan Dokkum et al.l (|2010T l . 

An alternative projection of Figure [7] is shown in Fig- 
ure [5] This figure shows the radius enclosing a given 
stellar ma ss as a function of red shift (analogous to Fig- 
ure 1 in iDiemand et al.l [2007) for galaxies at n c — 
1.4 x 10 -4 Mpc~ 3 . It is the equivalent of taking hori- 
zontal cuts in Figure [6l), and scaling by a stellar mass. 
The figure therefore depicts the evolution in the "onion 
layering" of stellar mass. The lines of constant mass ter- 
minate in high redshift bins where the given mass had 
not yet been assembled. Below z < 2, the radius enclos- 
ing a given mass remains roughly constant. This again 
highlights that new stellar mass growth below z < 2 oc- 
curs at larger radii and that the mass within the inner 
parts remains roughly unchanged. Between z ~ 3 and 
z ~ 2, the radius enclosing a given stellar mass appears 
to decrease by roughly a factor of ~ 2. 

5. DISCUSSION 
5.1. Inside-out Mass Growth 

In Section 0] we have shown that as a galaxy grows in 
stellar mass by a factor of ~ 3 from M ~ 5 x 10 10 M & 
at z = 2.75 to M ~ 1.5 x 10 11 M Q at z ~ 0, most of 
the new stellar mass that is added contributes towards 
mass growth at larger radii as one moves towards lower 
redshifts. The mass profiles presented in Figure [5] high- 
light this point as they show that most of the mass in 
the core is in place by z ~ 2. Below z ~ 2, Figure [7] 
shows how mass growth continues in the outer parts to 
the present epoch. As a consequence, the effective ra- 
dius of the galaxy grows from r ~ 2 kpc at z ~ 2 — 3 to 
r r~j 7 kpc at z ~ 0. Fitting the mass and effective radius 
evolution at z < 2, we find that r e oc M 2,0±0,3 , which is 
almos t exactly the relation found in Ivan Dokkum et all 
(pOlOh at z < 2. 

The high resolution of the HST imaging allowed us 
to probe the mass distribution of galaxies at z > 2 

at very small rad ii. This was not possible in the 

Ivan Dokkum et aLl (|2010f ) study as it relied on ground 
based imaging. With the higher resolution HST imag- 
ing, we find that the mass within r = 2 kpc, which is 
roughly the median effective radius for a galaxy selected 
at n c at z > 2, remains roughly constant at z < 2. Thus, 
there does not appear to be any substantial growth at 
z < 2 in the central part of the bulge that will charac- 
terize this galaxy at z ~ 0. From z ~ 3 to z ~ 2 there is 
evidence for mass growth in the inner part of the galaxy, 
as seen in Figures [7] and [H We note, however, that se- 
lecting galaxies at different values of n c can result in a 
more constant enclosing radius across redshift for a given 
mass. Larger samples at z > 2 may help in clarifying the 
potential mass buildup in the inner parts of galaxies at 
those early times. 

Dissipationless minor mergers are thought to play a 
significant role in the buildup of mass for QGs. In ad- 
dition, such processes are predicted to assemble mass at 
large radii, ther eby contributing towards the size growth 
of galaxies fe.g iBournaud et alj |2007t iNaab et aij|2009t 
iBezanson et al.l [20091) . For example, using a cosmologi- 
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cal hydrodynamical simulation, iNaab et al.l (|2009f ) sug- 
gest that az~0 galaxy with M ~ 1.5 x 10 11 M Q grew 
by a factor of ~ 3 after z ^ 3, primarily through dry mi- 
nor mergers. This mass growth is almost exactly what is 
found in this work for the same mass descendant at z ~ 0. 
Over this 11 Gyr timespan, the galaxy grew in size by 
a factor of ~ 3 — 4, also in agreement with our results. 
One caveat in this comparison is that the overall sizes 
of the simulated galaxy are a factor of 2 — 3 lower than 
what is found in this work. With more recent simula- 
tions, [Hiiz^eFlil] (|2012f ) sh ow that the relation r e oc M 2 , 
found in this work and in Ivan Dokkum et all (|2010D , is 
most easily reproduced by ~ 3 — 5 mergers with mass 
ratios of 1:5. In contrast, major mergers lead to a close 
to linear dependence of r e on M, which is not supported 
by the observations. It remains to be see n whether these 
mergers are actually obser ved (see, e.g.. iWilliams et all 
I20TQ INewman et al.1120121) . 

5.2. The Progenitors of Local Ellipticals at z ~ 3 

In the local universe, the most massive systems 
are generally bulge dominated elliptical ga laxies (e.g., 
Ivan der Wei et all 120091: IHolden et all 120121 ) with large 
ef fective radii (e.g., IShen et al.l 12003). In this work and 
in lvan Dokkum et al.l (|2010l ). we showed that the progen- 
itors of these local ellipticals at 1 < z < 2, are primarily 
compact QGs. As discussed in the previous section, these 
compact QGs are gener ally considered the c ores of ellip- 
ticals at z ~ (see also. lHopkins et al"]|2009l) . growing in 
mass and size potentially through dissipationless minor 
mergers to match the properties, namely sizes, of local 
ellipticals. 

With the HST imaging and the deep-IR UDS data, 
our number density selection allows us to the trace the 
properties of progenitors of local ellipticals to z ~ 3. At 
these redshifts, Figures [2] and [3] show that most of the 
progenitors are SFGs. The Sersic indices in Figure 0f> 
suggest these SFGs have exponential profiles, which is 
typically associated with the surface brightness profiles 
of disks. The Sersic index alone, however, is not definitive 
in defining the shapes of these progenitors at z ~ 3. In- 
stead, the most compelling evidence comes from the axis 
ratios shown in Figure 0b. At z ~ 0, galaxies selected at 
n c = 1.4 x 10~ 4 Mpc~ 3 have median observed axis ratios 
of b /a ~ 0.75. The imp lied intrinsic axis ratio is roughly 
2:3 (jHolden et a l. 2012), indicative of spheroidal systems. 
Meanwhile, at z ~ 3, Figure 0b shows that galaxies at n c 
have a median axis ratio of b/a ~ 0.52. This value for the 
axis ratio is very low considering that for a population of 
randomly oriented, infinitely thin disks, the median axis 
ratio would be b/a ~ 0.5. It is therefore likely that the 
progenitors of massive, quiescent, local elliptical galaxies 
at z ~ 3 are star forming thin disks. 

Many of the star forming disks at z ~ 3 disappear over 
the redshift range 1.5 < z < 3 and give way t o compact 
QGs, as seen in Figured iBarro et all (|2012f ) find that 
as the number density of extremely compact SFGs de- 
clines over 2 < z < 3, the number density of compact 
QGs rises in a complementary way, suggesting an evo- 
lutionary link. A scenario in which star forming disks 
simply quench through some mechanism and passively 
evol ve to become compact Q Gs may not be feasible given 
that [Whitakcr et all (|2012l) find that the most recently 
quenched galaxies at any time have compact sizes that 



are smaller than those of SFGs. 

6. SUMMARY 

We used HST imaging to study the structural proper- 
ties of galaxies selected at a constant cumulative num- 
ber density of n c = 1.4 x 10 -4 Mpc -3 at redshifts of 
0.25 < z < 3. This selection allowed us to trace the evo- 
lution of galaxies with stellar mass M = 5 x 10 10 M@ 
at z = 2.75, as they grew by a factor of ~ 3 to become 
~ 2M* galaxies in the local universe. This work build s 
on the previous analysis by Ivan Dokkum et all (2010), 
who also selected galaxies at a constant number density. 
Here, we employ high resolution HST imaging and ex- 
tend the analysis to z ~ 3. In contrast to mass-selected 
studies, our selection at a constant cumulative number 
density allows for a more straightforward evolutionary 
link between progenitors and descendants at different 
redshifts. At 1< z < 3, we employed CANDELS WFC3 
J125 and -H160 imaging in the UDS, while at 0.25 < z < 1 
we used wide-field ACS Igu imaging in COSMOS to fit 
single-component Sersic profiles at a common rest-frame 
wavelength. The resulting Sersic indices, effective radii, 
and axis ratios were used to aid in our analysis. The uni- 
form datasets and analysis methods carried out in this 
work serve to minimize systematics. 

Our main conclusions are the following: 

1 . The typical galaxy at the selected value of n c has 
grown in effective radius by a factor of ~ 3 — 4, 
mostly since z ~ 2 (Figure 0]) . 

2. The evolution in the stellar mass surface density 
profiles of galaxies selected at n c indicate that most 
of the mass in the central regions was in place by 
z ~ 2, while almost all of the new mass growth 
took place in the outer parts (Figures [6j [3 and 
[5]). This inside-out mass growth is responsible for 
the increase in size and Sersic index towards low 
redshift. 

3. At z < 2, we find that as the stellar mass builds up, 
the effective rad ius grows as r e oc M 2 °, in excellent 
agreement with Ivan Dokkum et"aH (|2010D . Recent 
simulations show that such a dependence is con- 
sistent with mergers with mass ratios o f 1:5 being 
respo nsible for much of the size growth (|Hilz et al.l 
12011) . 

4. At z ~ 3, the rest-frame UVJ colors, Sersic in- 
dices, and axis ratios indicate that the progenitors 
of present day massive galaxies were star forming 
disks with r e w 2 kpc and a third of the z ~ stel- 
lar mass. At 1.5 < z < 2, these galaxies doubled in 
stellar mass and were mostly compact QGs. These 
galaxies evolved further into the ~ 2M* galaxies in 
the local universe that are known to be quiescent, 
bulge dominated, elliptical galaxies with large ef- 
fective radii. 
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APPENDIX 

THE STRUCTURAL PROPERTIES OF QUIESCENT AND STAR FORMING GALAXIES FOR STELLAR MASS LIMITED 

SAMPLES 

The primary purpose of the paper was to examine the structural properties of galaxies selected at a constant 
cumulative number density. In doing so, we were able to the trace evolution of ~ 2M* galaxies from their progenitors 
at z ~ 3. Here we show results for galaxies selected above a constant mass limit, a more commonly used selection 
that provides an alternative view to the number density selection. We distinguish quiescent galaxies (QGs) and star 
forming galaxies (SFGs) with a L/FJ-selection, as these two subpopulations are known to have differing properties. 
As the mass selection allows for a larger sample, we can confirm some of the general findings in the main part of the 
paper where the sample was much smaller due to the number density selection. 

Classification of QGs and SFGs with UVJ Selection 

Figure |H] shows UVJ diagrams for galaxies in different redshift bins at0.25<z<3 abo ve a stellar mass o f 
M > 10 10 5 Mq. Above this mass limit, the sample is complete for QGs and SFGs to z ~ 2.5 (Quadri ct al.|[20T2T) . 
thus the highest redshift bin at 2.5 < z < 3 likely exhibits some incompleteness for QGs. The subset of objects 
with measured structural parameters from the HST imaging are color-coded according to their Sersic index. For the 
COSMOS sample at 0.25 < z < 1, each color-color bin is color-coded according to the median Sersic index for galaxies 
in that bin and the symbol size reflects the size of the sample within the bin (on a logarithmic scale). The boundary 
distinguishing QGs from SFGs is shown for each redshift bin. 

For the redshift range studied here, Figure [S] shows the buildup in the proportion of QGs, down to a fixed mass limit, 
with time. The figure also shows that QGs generally have higher Sersic indices relative to SFGs at a given redshift. 
This is also the case for optically red galaxies (e.g., U — V > 1.5), where C/T^J-selected QGs have higher Sersic indices 




0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 

V-J (AB mag) V-J (AB mog) V-J (AB mog) 

Figure 9. Rest-frame U - V vs. V - J for galaxies with mass, M > 10 10 5 Mq, in COSMOS and the UDS. In COSMOS (0.25 < z < 1), 
each color-color bin is color-coded by the median Sersic index and the symbol size is indicative of the number of galaxies in the bin (weighted 
on a logarithmic scale). Bins with less than two galaxies are not color-coded. In the UDS (1 < z < 3), the gray points indicate the parent 
sample of galaxies while objects color-coded by Sersic index indicate the main subsamplc with HST imaging. The sample size of galaxies 
with measured structural parameters is indicated in the bottom right of each panel. The solid black line shows the division between QGs 
and SFGs. At a given redshift, UV J-selected QGs generally have higher Sersic indices compared to SFGs. 
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relative to reddened SFGs. The UVJ selection therefore works efficiently to distinguish galaxies to z ~ 3 based not 
only on their SFHs, but also their structural properties. This was suggested in the analysis at constant number density, 
but is confirmed here with a much larger sample. The fact that QGs and SFGs have different structural properties 
at the highest redshifts provides additional confidence in the structural parameters obtained from GALFIT: high and 
low Sersic index galaxies are found where we think they should lie, in the quiescent and star forming regions of UVJ 
color space, respectively. 

Sersic Indices 

Figure [TO] shows the median Sersic index, n, of QGs and SFGs at different redshifts. Error bars (±lcr) were computed 
by bootstrapping the sample in different redshift bins. Note that the random errors at z < 1 are generally smaller 
than the size of the data point given the large sample size in COSMOS. QGs (red circles) show mild evolution in their 
Sersic index over 0.25 < z < 3, increasing slightly from n s» 3 to n s» 4 at low redshift. A fit of the form n cx (1 + z) a 
indicates that a = —0.50 ±0.18. We ignore the z ~ 2.75 QG data point for this fit as incompleteness likely plays a role 
for QGs at the highest redshifts. In carrying out these fits, we also add a 10% systematic uncertainty in quadrature 
to the random errors to account for potential differences between the COSMOS and UDS samples. The median Sersic 
index for SFGs (blue circles) increases from n«lat2<z<3ton«2at low redshift. A fit of the form n cx (1 + z) a 
indicates that a = — 0.64± 0.16. The Sersic indices of both QGs and SFGs with mass M > 10 10 5 M are therefore 
decreasing towards higher redshifts. 

A potential technical explanation for the declining Sersic indices of QGs with redshift is the presence of SFGs which 
scatter into the QG bin. The rest-frame colors are more uncertain at higher redshifts and there are relatively more 
SFGs, which generally have lower Sersic indices. We test this explanation by comparing the SSFRs of QGs with high 
Sersic indices (n > 2.5) with those that have low Sersic indices (n < 2.5). Using the SSFRs from the SED fit we find 
that for the different redshift bins above z > 1.5, the QGs with low and high Sersic indices generally have median 
SSFRs within the uncertainties. We note that for the highest redshift bin, 2.5 < z < 3, the difference in SSFR is 
somewhat larger. However, incompl eteness for QGs l ikely plays a role for this redshift bin. It is worth noting that for 
a similar mass limit as in this work. [Newman et al.l (|2012t ) also find that the Sersic indices of QGs decline to z = 2.5 
with QGs selected based on an SSFR limit and a lack of detection in MIPS 24 fim imaging. Likewise for SFGs at low 
redshift, the more numerous QG population could scatter into the SFG selection, leading to an elevated median Sersic 
index for SFGs. The smaller uncertainties for the rest- frame colors make this scenario less likely. We confirm that the 
typical SSFRs of SFGs at 0.25 < z < 1 with low and high Sersic indices are the same within the uncertainties. We 
therefore conclude that the decline with redshift in the Sersic index for QGs and SFGs with mass M > 10 10 5 M is 
not a consequence of galaxies of either type scattering into the other bin. 

Again, we note that at the highest redshifts, the Sersic indices obtained from GALFIT are generally different for 
QGs and SFGs, as is the case at lower redshifts where the signal-to-noise of the measurements are more robust. In 
addition to our simulations discussed earlier, this point affirms that our structural parameters are reasonable at high 
redshift. 
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Figure 10. Sersic index vs. redshift for QGs (red) and SFGs (blue) with stellar mass M > 10 10 - 5 Mq. The la error bars are computed 
from bootstrapping. Dashed lines represent a fit to the data of the form n oc (1 + z) a , where a = — 0.50 ± 0.18 and —0.64 ±0.16 for QGs and 
SFGs respectively. For QGs, the highest redshift data point at 2.5 < z < 3 (open circle) is not included in the fit due to incompleteness. 
Because of potential systematic differences in the COSMOS (0.25 < z < 1) and UDS (1 < z < 3) samples, the fits were carried out with 
a systematic uncertainty of 10% added in quadrature to the la errors shown. The Sersic indices of both QGs and SFGs decrease towards 
higher redshift. 
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Figure 11. (a) Effective radius vs. redshift for QGs (red) and SFGs (blue) with stellar mass M > 10 Mq. The la error bars are 
computed from bootstrapping. For a galaxy at a given redshift, the HST bandpass nearest to rest-frame ~ 5200 A was selected for the size 
measurement. Dashed lines represent a fit to the data of the form r e <x (1 + z) a , where a = —1.16 ± 0.20 and —0.63 ± 0.13 for QGs and 
SFGs respectively. For QGs, the highest redshift data point at 2.5 < z < 3 (open circle) is not included in the fit due to incompleteness. 
The sizes of both QGs and SFGs decrease towards higher redshift. (b) Mass- normalized effective radius vs. redshift. The effective radius 
of each galaxy has been scaled to the value for that of a 10 11 M p, galaxy (Mu = Af/10 11 Mq) assuming that r e oc M@ at all redshifts, 
where /3 = 0.56 for QGs and /3 = 0.3 for SFGs (Shcn ct al. 2003). The redshift evolution for the mass-normalized radius is characterized 
by a = -1.30 ± 0.20 and a = -0.80 ± 0.13 for QGs and SFGs, respectively. 

Sizes 

Figure [Tib shows the median effective radius (r e ) of QGs and SFGs at different redshifts. QGs above M > 10 105 Mq 
show substantial evolution in their sizes over 0.25 < z < 3, increasing by roughly a factor of ~ 3 from r e ~ 1 kpc to 
r e ~ 3 kpc. A fit of the form r e oc (1 + z) a indicates that a = — 1.16 ± 0.20 for QGs. The QG data point at 2.5 < z < 3 
was again not included in the fit due to incompleteness. Meanwhile, SFGs increase in size by a factor of ~ 2 over the 
same redshift range from r e ~ 2.5 kpc to r e ~ 4.5 kpc. We find a = —0.63 ± 0.13 for SFGs. We note that at high 
redshifts, the median sizes for QGs and SFGs are always larger than the measurement limits imposed by the PSF of 
the WFC3 J125 (FWHM/2 - 0'.'06) and H 160 (FWHM/2 - 0'.'09) imaging. 

Figure [Tib shows the median mass- normalized effective radius versus redshift for QGs (red) and SFGs (blue). The 
effective radii were normalized to a stellar mass of 10 11 Mq assuming a size-mass relation r e oc M@ , where j3 — 0.56 
for QGs and (3 = 0.3 for SFGs (|Shen et al.ll200l . The mass-normalized radii therefore represent the size that galaxies 

would have if they lie on the given size-mass relation with a stellar mass of 10 11 Mq. A fit of the form r e jM^ x oc (l+z) a , 
where M n = M/10 11 Mq, indicates that a = -1.30 ± 0.20 for QGs and a = -0.80 ± 0.13 for SFGs. 

For our mass-limited sample, the median size of QGs is generally smaller than that of SFGs at a given redshift 
and the sizes of both QGs and SFGs decrease towards higher redshift as also found by several other authors (e.g ., 
iDaddi et al.lf2005t iTofTeTaLl [20071: ICimatti et al.ll200l IPranx et al.ll200l iWilliams et ai.|[2010t iNewman et al.lf2012f ). 

Axis Ratios 

The distribution of axis ratios provide additional insight into the shapes of QGs and SFGs. Figure [12] shows the 
median axis ratio of QGs and SFGs with mass M > 10 5 Mq at different redshifts. The dashed lines are fits of the 
form b/a oc (1 + z) a . For QGs, a fit to the data at 0.25 < z < 2.5 indicates a — —0.01 ± 0.17, consistent with no 
change over this redshift range in the media n b/a. The med i an ax is ratio for all QGs at 0.25 < z < 2.5 is b/a ~ 0.69, 
which is very close to the value computed by INewman et al.l ()2012[ ) for a similar mass limit (after accounting for IMF 
differences) and redshift range. Incompleteness impacts the QG data point at 2.5 < z < 3 and is the refore ignored 
in the fit above. Also shown in Figure IT2l is the median b/a value at z ~ from SDSS determined bv iHolden et al.1 
(2012) for galaxies with stellar mass 10 10 5 < M/Mq < 10 11 . IHolden et all (|2012j ) find no change in the distribution 
of axis ratios of QGs from z ~ to z ~ 0.7, consistent with our results at those redshifts. For a sample of 14 QGs 
at 1.5 < z < 2.5 with stellar masses M > 10 10 ' 8 Mq, Ivan der Wei et all ()2011h computed a median axis ratio of 
b/a = 0.67. For the same redshift ra nge and mass cut we m easure a median axis ratio of b/a — 0.73 ± 0.06 for 46 
QGs, which is consistent with that of Ivan der Wei et al.l (|201 If ) . They suggest that QGs may be more disk dominated 
at high redshift based on the proportion of QGs with low axis ratios. A more detailed analysis of the distribution of 
axis ratios with larger samples may lead to a more definitive answer. 

The median axis ratio of SFGs also does not change significantly over 0.25 < z < 3. We find a = 0.01 ± 0.13, 
consistent with no evolution. The median value for the axis ratio of SFGs over 0.25 < z < 3 is b/a ~ 0.58. The 
deviation at 1.5 < z < 2 is within ~ 1.8<7 of this value. Interestingly, for a population of infinitely thin disks, the 
median observed axis ratio would be b/a ~ 0.5, not far from the value found here. 

The mass selection affords a larger sample at 2 < z < 3 than was possible with the number density selection. We find 
the general trends found with the latter selection still hold with this much larger mass limited sample. In particular, 
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Figure 12. Median axis ratio, b/a, versus redshift for QGs (red) and SFGs (blue). The 1<t error bars are computed from bootstrapping 
each subsam p le. T he filled circle at z ~ represents the median b/a value of QGs with mass 10.5 < log M/Mq < 11 computed by 
IHolden et al.f 120121) . Dashed lines indicate fits to our data of the form b/a tx (1 + z) a . The QG data point at 2.5 < z < 3 (open circle) 
was not included in the fit due to incompleteness in that bin. For both QGs and SFGs, a is consistent with no evolution in the median 
axis ratio. 

the axis ratios of SFGs, which dominate both the mass limited sample and the number density selected sample at 
z > 2, display median axis ratios consistent with randomly oriented disks. 



Stellar Mass Variations 

The implemented stellar mass limit, M > 10 10 5 M Q , to first order limits variations in the distribution of stellar 
masses across different redshift bins. This works to minimize residual correlations between various parameters and 
stellar mass. For QGs, the range of median stellar masses for the different redshift bins is only 0.07 dex. For SFGs, 
the range in median stellar masses is somewhat higher at 0.18 dex, with the higher redshift bins (z > 1.5) having 
more massive SFGs above the mass limit. We do not expect this difference for SFGs to have a significant impact on 
our conclusions. In fact, it strengthens many of our key points. For example, if higher mass SFGs have higher Sersic 
indices, then the Sersic indices for SFGs at 1.5 < z < 3 should be even lower than those at low redshift after correcting 
to the same median mass. 



